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At doping below 6% the bilayer cuprate YBa2Cu306+i/ is a coUinear antiferromagnet. Independent 
of doping the value of the staggered magnetization at zero temperature is about 0.6^s. This is the 
maximum value of the magnetization allowed by quantum fluctuations of localized spins. In this 
low doping regime the compound is a normal conductor with a finite resistivity at zero temperature. 
These experimental observations create a unique opportunity for theory to perform a controlled 
calculation of the electron spectral function. In the present work we perform this calculation within 
the framework of the extended t — J model. As one expects the Fermi surface consists of small 
hole pockets centered at (±7r/2, ±7r/2). The electron spectral function is very strongly anisotropic 
with maximum of intensity located at the inner parts of the pockets and with very small intensity 
at the outer parts. We also found that the antiferromagnetic correlations act against the bilayer 
bonding-antibonding splitting destroying it. The bilayer Fermi surface splitting is practically zero. 

PACS numbers: 74.72.Gh, 75.10.Jm, 75.50.Ee 



I. INTRODUCTION 

One of the central issues in the physics of cuprates is 
the evolution of Fermi surface upon doping. Despite a 
consensus on the large Fermi surface in the overdoped 
side, the shape and nature of the Fermi surface in under- 
doped cuprates remains an unresolved issue. Cuprates 
are doped Mott insulators. Theoretically there is no 
doubt that minima of the dispersion of a single hole 
injected in the Mott insulator are at the nodal points 
(±7r/2, ±7r/2). This conclusion is supported by the 
Angle-Resolved Photoemission Spectroscopy (ARPES) 
data from undoped cuprates.— This implies that at a 
sufficiently small doping holes must go to the hole pock- 
ets. This is the small Fermi surface situation. In real 
compounds disorder complicates the situation and signif- 
icantly masks the generic physics. For example the proto- 
typical cuprate La2-xSra;Cu04 is so disordered that holes 
are strongly localized at doping x < 0.1. In the case of 
strong localization the notion of the Fermi surface is am- 
biguous. ARPES data from underdoped La2-i^Sr^Cu04 
reveal Fermi arcsii^ that are not consistent either with 
the small (hole pockets) nor with the large Fermi sur- 
face. On the other hand the very recent ARPES data 
from underdoped Bi2Sr2CaCu208+5 indicate small hole 
pockcts;2, 

YBa2Cu306+y (YBCO) is probably the least disor- 
dered cuprate in the low doping regime. In this paper 
we denote doping x to be the hole concentration per unit 
cell of the Cu02 layer. YBCO is not superconducting 
below X ^ 0.06, where the compound remains a normal 
conductor with delocalized holes. The zero temperature 
resistivity remains finite)^ apart of a very weak logarith- 
mic temperature dependence^i^ expected for a weak dis- 
order. The heat conductivity also indicates delocalization 
of holes ii We emphasize that this is very different from 
La2-2;Sr2;Cu04 where holes are localized and hence the 



compound is the Anderson insulatoriSiS at x ^ O.lJ^ Ul- 
timately, at the very low doping, x ^ 0.01, the disorder 
wins even in YBCO and it also becomes the Anderson 
insulator)^ It is helpful to have in mind an approximate 
empiric formula^iii 

X « 0.35(y- 6.20) (1) 

to relate the doping level x and the oxygen content y in 
underdoped YBa2Cu306+y at x < 0.12. The homogene- 
ity of YBCO is the reason why Magnetic Quantum Oscil- 
lations (MQO) were observed in this compoundii2r^ The 
oscillations clearly indicate small Fermi pockets, while 
strictly speaking nature of the pockets, including the sign 
of the charge of the fermion, experimentally remains a 
controversial issue. 

The collinear antiferromagnetism (AF) in YBCO is 
preserved up to the doping level x ~ 0.06. Moreover, 
the zero temperature staggered magnetization ^ « 0.6^b 
is practically doping- independent, having the same value 
as in the parent Mott insulator This is the maximum 
value of magnetization allowed by quantum fluctuations 
of localized spins. The doping behavior of the staggered 
magnetization in YBCO is very different from that in 
La2-2:Sr2;Cu04 where the staggered magnetization de- 
cays dramatically with doping. A special mechanism has 
been proposed^^ to explain antiferromagnetism in YBCO 
at a; ^ 0.06. For purposes of the present work details of 
the mechanism are not important. The only important 
point is that the antiferromagnetism is independent of 
doping. This direct experimental observation in combi- 
nation with simple metallic behaviour (finite resistivity 
at zero temperature that is another direct experimen- 
tal observation) gives a unique opportunity to perform 
a controlled and fully reliable theoretical calculation of 
the electron spectral function at x ^ 0.06. In the present 
work we perform this calculation. 

Our analysis is based on the t — t' — t" — J model and 
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employs the self-consistent Born approximation (SCB A) . 
SCBA has been widely applied to study a single hole 
dressed by spin fluctuations.—"— On the other hand ap- 
plication of the method at finite doping has been very 
limite d^^'^° because a usual justification of the method 
requires a long range AF order, and the common wisdom 
is that even a tiny doping destroys the order. AF ordered 
YBCO docs not comply with the common wisdom and 
provides a unique opportunity to address spin fluctua- 
tions up to very high accuracy. The AF order implies 
that the single loop vertex correction is forbidden (the 
"Migdal theorem"), so SCBA is exact up to double loop 
corrections 3^ 

Another important property of the quasiparticle dis- 
persion is the layer bonding-antibonding splitting in the 
bilayer cuprate. A density functional theory calculation 
shows a more or less constant splitting through out the 
whole Brillouin zonci^ ARPES in overdoped YBCO con- 
firms the splittingi^i^ On the other hand, in the under- 
doped regime ARPES measurements^^ indicate no such 
splitting. To explain this in the present paper we show 
that antiferromagnetic correlations between layers in the 
bilayer system diminish the splitting in spite of the strong 
chemical tendency towards the splitting. The similar 
result for undoped YBCO was previously obtained in 
Ref. [13. The antiferromagnetic correlations are due to 
the antiferromegnetic exchange J± ~ lOmeV between 
the layers.— It is worth noting that the small hole pock- 
ets are essential for the suppression of the bilayer splitting 
in the underdoped regime. 

In the present work we consider the clean limit without 
any disorder. There are recent ARPES experiments on 
underdoped YBCO^i^ where the mechanism of doping 
is related to depositing of potassium atoms on the sur- 
face. The potassium ions give rise to a random potential 
for mobile holes. The degree of disorder in these exper- 
iments remains an open issue when compared with the 
present results. 

The structure of the paper is the following. In Section 
II, we introduce the Hamiltonian and calculate vertexes 
for many-body diagrammatic technique. In Section III 
the method of SCBA at finite doping is discussed. The 
method implies the spin-charge separation. The spin- 
charge recombination amplitude and the ARPES spectral 
function is calculated in Section IV. Section V summa- 
rizes our results. 



II. THE EXTENDED t J MODEL FOR THE 
BILAYER YBCO 



We simulate lightly doped YBCO by using the double 
layer t—t'^t" — J model with constant interlayer coupling 
t± and J± 



La 



— Hfj" + Hj^jj^ + Ht. 



(2) 



where t/t'/t" is the nearest /next- nearest/next-next- 
nearest neighbor in-plane hopping, respectively, and m = 
{1,2} is the plane index. Since J± ~ lOmeV > 0, the 
magnetic ordering along c— axis is also AF. We denote 
the projected coordinate in each plane as 

i G t sublattice plane 1=1 sublattice plane 2 
j G i sublattice plane 1 = ^ sublattice plane 2 
l^ihj} (3) 

Through out the article we set energy unit as J = 
ISOmeV— >■ 1, hence J± w 0.08. Schematics of the model 
is shown in Fig. [T] 




plane 1 



plane 2 



J 



rn,{ij)^a rn^{ij)' ,a 



FIG. 1: (color online) Schematics of double layer Hamiltonian, 
Eq. and definition of coordinate, Eq. ([3]). 



We use the following values of the intralaycr hopping 
parameters t = 3.1, t' = -0.5, t" = 0.4. While the 
value of t corresponds to that obtained in the density 
functional theory calculation^ values of t' = —0.5 and 
t" = 0.4 are somewhat different. Ref. [Sll gives for the 
optimally doped YBCO the following values, t' « 0.8, 
t" w —0.7. Based on our results that we compare with 
ARPES data we believe that the values of t' = —0.5 and 
t" = 0.4 accepted in the present work are more suitable 
for underdoped YBCO. However, in the end the differ- 
ence between these two sets of parameters is not qualita- 
tively important. We take a constant interlayer tunneling 
t± based again on the first principle calculation)^ which 
shows a practically constant splitting between bonding 
and antibonding bands with the value of the splitting 
corresponding to t± sa lOQmeV 1. ARPES data from 
overdoped YBCO^^i^^ support this splitting. The esti- 
mate « 1 is supported also by the ratio of the su- 
perexchange parameters, 0.08 = J±/J = t\/t^- 

Following the standard SCBA philosophy the Hamilto- 
nian ^ is grouped into three sectors, Hf^t", Hj j^, and 
Ht^t± , which correspond to bare hole, bare magnon, and 
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hole-magnon interaction, respectively. We first discuss 
the bare magnon sector Hj j^. The Holstein-Primakoff 
bosons are defined according to sublattices in each layer, 
as defined in Eq. ([3]). In plane 1 

_ 1 t 
— "1,1 



(4) 



and in plane 2 



Sj 2 = ^1,2 



^3,2 - 2 " aj,2aj,2 



t 

The Fourier transform is 



(5) 



(6) 



where summation over q is restricted inside Magnetic 
Brillouin Zone (MBZ). By introducing parity + and — 
bases with respect to interchange of two planes 



(7) 



and Bogoliubov transformation 



^q,± = "q,±"q,± + ^'q,±/3lq,± 



^-q,± = «q.±aq,± + Mq,±/3-q,± (8) 

one can diagonalizc the Hamiltonian^ii^ 



q 



7q = ^ (cosq.^ + cos^j^) 



(9) 



where we denote a_L = J±/ J- Here Wq^+ mode is gapless 
at (0,0) and gapped at (7r,7r), while Wq._ mode is the 
opposite. This is consistent with the fact that optical 
mode is frequently referred to even parity, and acoustic 
mode to odd parity in inelastic neutron scattering experi- 
ments which practically measure magnon dispersion near 
(7r,7r)i^i^ The Bogoliubov coefficients are 



Uq,± 




Wq,± = -sign(7q ± , , 

4 y Wq,± 

The magnon Green's function is defined as 
Jo 

= / (T/3q,±(t)/3t ±(0))e-*dt 



(10) 



(11) 



In the present work we do not consider renormalization 
of magnon dispersion due to interaction with holes, but 
simply adopt the magnon sector of undopcd Mott insula- 
tor. This assumption is justified because /xSR measure- 
ments indicate almost unrenormalized staggered magne- 
tization up to doping level x ~ 0.06.— The magnon dis- 
persion close to the Goldstone point q = [equivalent 
to q = (±7r, ±7r)] is somewhat changed under dopingi^ 
However, due to the Adler's theorem, magnons close to 
the Goldstone point practically do not influence the hole 
dispersion, see discussion below. The main contribution 
to the hole dispersion comes from magnons that are far 
away from the Goldstone point. For this regime, res- 
onant inelastic X-ray scattering(RIXS) clearly demon- 
strates that magnons are almost independent of doping.— 
Theoretical analysis of the RIXS experimental data is an 
important problem that will be considered elsewhere.— 
To address the bare hole dispersion from the -fft'.t" 
term, we first define hole operators according to the co- 
ordinate in Eq. ([3]) 



dk.i, 



N{l/2 + m) ^^^-^'^ 



Ed.. 



N{l/2 + m) 



E4.- 
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where m = 
states are 



iV(l/2 + m) 



iV(l/2 + m) 



E4 



(12) 



1(5",) I w 0.3, see Ref. The fixed parity 



la) 



(b) 



(5, k V, k+q 



+ + .'' ^\ + + 



+ 

+ 



(13) 



(c) 



The bare hole dispersion is then 



ej^ = 4t' cos kx cos ky + 2t" (cos 2kx + cos 2A:y) 

(14) 



Since site i on plane 1 and site i on plane 2 belong to dif- 
ferent sublatticcs, the interlayer hopping does not enter 
the bare dispersion. This is the reason why the AF cor- 
relations suppress the interlayer hopping. Similar to the 
in-plane nearest-neighbor hopping, the interlayer hop- 
ping contains the spin-flip process that contributes to 
the hole-magnon interaction. The corresponding vertex 
therefore contains contributions from both t and t±, as 
calculated in Appendix A. The Hamiltonian -fft,*^ and 
the corresponding vertex read 

k,q,7(5i/ 

4 



(15) 



+4+q.,.,tC'k,54/5q,7 ) + h-C. , 
5k,q,75 = —j= (7k'"q,7 -I- 7k+qfq,7) 

V J V 



■ 7Wq,7) 



We denote the parity index by 7, (5, = ±1. This is the 
parity of the annihilated magnon, the annihilated hole, 
and the created hole, as shown in Fig.[2][a). Conservation 
of parity implies v = j6. The vertex psp is zero at 5 = 0, 
this is a consequence of Adler's theorem, and this is why 
the contribution of magnons with small momenta to the 
self-energy is negligible. 



III. SELF CONSISTENT BORN 
APPROXIMATION 

The SCBA is based on summation of diagrams having 
the highest possible power of the hopping parameter t 
at a given number of loopsi^"— This implies that the 
method is justified at t/J > 1. A very important point 
is that due to the spin structure of the theory (conser- 
vation of Sz) there is no single loop correction to the 
hole-magnon vertex shown in Fig. ^a). Therefore, the 
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FIG. 2: (color online) Feynman diagrams for (a) the hole- 
magnon vertex defined in Eq. (|15|l . (b) Dyson's equation 
Eq. (|19p for Gd± , and (c) schematics of advanced and retarded 
part of each self-energy diagram, as calculated in Eq. (|Bip . 



diagrams having the highest possible power of t are the 
rainbow diagrams. Thus, SCBA is summation of rain- 
bow diagrams. Again, the method is justified due to (i) 
absence of the single loop vertex correction, (ii) due to 
t/J> 1. 

Compared to the single hole in a single layer 
case?^"— the present calculation has two complica- 
tions: the double layer and the finite doping. Impor- 
tantly, since we consider the coUinear AF state, the both 
points (i) and (ii) justifying the method are still valid in 
spite of the complications. We stress that although the 
finite doping case has been discussed before^^i^ it is not 
until current calculation that SCBA is rigorously formu- 
lated for a real material that AF order persists at finite 
doping. 

At finite doping we have to use the Feynman Green's 
function of the hole 



/oo 
-00 



(0))e'^*di (16) 



instead of the retarded Green's function in the undoped 
casei^"— The Green's function has the parity index 
± reflecting the double layer structure, and the pseu- 
dospin index a =ti reflecting the AF structure. Since 
the up and down pseudospins are degenerate, we omit 
the psedospin index a for the rest of the article, although 
one should keep in mind that vertexes in Eq. p5|) always 
flip the pseudospin. In the calculation of self-energy, we 
adopt the spectral representation 



Gd±(e,k) ^ / dx ^—77+ / 

Jo e-x + iO 



X — iO 



(17) 



The chemical potential is set equal to zero. The techni- 
cal advantage of using spectral representation is that it 
ensures the causality of Dyson's equation in any order, 
as we calculate the self-energy in terms of the spectral 
functions. This is important because Dyson's equation 
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typically converges at about 10 '-^ 20tli order, therefore 
causality must be ensured at each order. The doping is 

2a; = 2^^/ dxB±{x,k) . (18) 
k ± 

The summation over k is limited inside the MBZ. The 
coefficient 2 in the left-hand-side of ((T8)) is due to the 
bilayer, the coefficient 2 in the right-hand-side is due to 
pseudospin. 

Dyson's equations shown schematically in Fig. [2{b) 
read 

Gd+ (e, k) = [e - el - S++ (e, k) - (e, k) + iO] 

Gd-{e, k)=[e- el - S]+_(e, k) - k) + zO] 

(19) 

The subscripts in the self-energy S±,± indicate parities 
of the intermediate hole and magnon as it is shown in 
Fig.mjb). Note that the considered theory does not pos- 
sess the usual cross-leg-symmetry of vertexes. Therefore 
the usual Feynman technique is not sufficient)^ and one 
has to adopt the Goldstone-Brueckner technique that ex- 
plicitly separates forward in time and backward in time 
diagrams, as shown in[2l^c). Explicit expressions for the 
self-energy are derived in Appendix B and presented in 
Eq. (IbTI) . 

We now discuss symmetry properties of Green's func- 
tions imposed by definitions of operators ([7]) and ([T^. 
Due to the chequerboard AF order, hole and magnon op- 
erators change parity under translation by the AF wave 
vector Q = (tt, tt), k — > k + Q. 

aq,+ = "q,- 

/3q,+ = -/3q.- . (20) 

Hence the Green's functions satisfy the following symme- 
try relations 

Gd+(e,k + Q) = Gd_(e,k) , 

D+{cu,q + Q)=D^{io,q) . (21) 

As a result, the hole dispersion of parity + and — swap at 
the MBZ boundary, as previously reported for the single 
hole casci^ This implies that at the MBZ boundary, the 
magnon dispersion of either parity are degenerate, and 
the hole dispersion of either parity are also degenerate, 
as addressed below. 

Iterative numerical solution of Eqs. ([T^ . (|Bip requires 
more computational power compared to solution of sim- 
ilar equations for undoped single layer casci^"— Never- 
theless the solution is straightforward. We solve Eq. ^T9\\ 
inside MBZ on a 64 x 64 cluster, with about 1700 fre- 
quency points (the frequency grid is Aw ^ 0.015). We 
present plots of the hole spectral function Sd defined as 
A{e, fc) if e > and as B(e, fc) if e < 0, see Eq. p7| . 

Sd±{e,k) = A±ie,k)+B±ie,k). (22) 




FIG. 3: Spectral functions Sd+ (a) and Sd- (b) for different 
values of k along the nodal direction from (0,0) to (7r,7r). 
There is an offset 0.2 between subsequent values of k. The 
doping level is a; = 0.059. 



Figs.El^a) and (b) display the hole spectral functions (off- 
set 0.2) of positive and negative parities for doping level 
X = 0.059. This doping is close to the highest doping level 
at which AF order persists while superconductivity has 
yet taken placci^ To show more details the spectral func- 
tion 5'd+ is plotted inFig.Hfork = (0,0), k= (7r/2,7r/2) 
and k = (7r,0). The spectral functions show a striking 
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FIG. 4: Spectral functions Sd+ for k = (0,0), k = {tt/2,tt/2) 
and k = (tt, 0). The doping level is a; = 0.059. 



similarity to the hole spectral function in an undoped in- 
sulator: there is a well defined quasiparticle peak at low 
energy, and a large incoherent part at higher energy. In 
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addition, a small but significant incoherent part is ob- 
served at energy below the quasiparticle peak, due to 
the electron plus multiple magnon configurations in the 
emission channel. Along nodal direction, the quasiparti- 
cle peak is most pronounced near (7r/2,7r/2), and grad- 
ually decreases as approaching (0, 0) and (tt, tt) where it 
is dissolved in the hole-magnon continuum. 

The position of the quasiparticle peak determines the 
quasiparticle dispersion e^. The dispersion is shown in 
Fig.[5ja) along a certain line in the Brillouin Zone (BZ). 
Formation of small Fermi surface with hole pockets at 
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FIG. 5: (a) The quasiparticle dispersion, and (b) the quasi- 
particle residue Zd along a certain line in BZ. The doping 
level is a; = 0.059. 



k = (7r/2,7r/2) is evident. Eq. (PO]) . imposes a nontrivial 
constraint on the dispersion: From k = (0,0) to (7r,7r), 
one sees that the two parity bands are symmetric but 
swap at the MBZ boundary k = (7r/2,7r/2). Along the 
MBZ boundary, as shown in the section from k = (tt, 0) 
to (7r/2,7r/2), the two bands are degenerate. At the 
Fermi surface ek+ differs from ek_ by only « 25meV. 
We remind that for the value t± = J = 130meV used in 
the calculation, one expects the splitting to be 260mcV. 
Therefore AF correlations strongly suppress the bilayer 
splitting. 

In the section of (7r,7r) to (7r,0) in Fig. [Sja), one sees 
a local dispersion minimum near (7r,7r/2). However, in 
the (7r/2,7r/2) to (7r,7r/2) section, ones sees that (7r,7r/2) 
is a local maximum. Hence, there is a saddle point near 
(tt, 7r/2), or cquivalcntly near (0, 7r/2). (7r/2, 0). The sad- 
dle point gives a very large contribution to the density of 
states, iV±(e) = J [A±{e,k) + B±{e,k)] d:^k/{2nf . Due 



to the saddle point the electron response function is very 
strongly peaked at q « (7r/2,0), ui « lOOmeV. Inter- 
estingly, these parameters are very close to where the 
anomaly in the breathing phonon mode is observedi^ 
Although further investigation is necessary to clarify this 
point, we anticipate that the saddle point plays a crucial 
role in the anomaly. 

The quasiparticle residue Zd]i± is shown in Fig. [Sfb) 
along the same line in BZ as dispersion. Along the nodal 
direction, the quasiparticle residue Ziik.± is maximum at 
the bottom of the pocket, and decreases smoothly as 
moving away from the pocket. This smooth dependence 
is very similar to the residue in undoped parent com- 
poundi^ Since the pocket is rather small, the residue in- 
side the pocket can be well approximated by a constant 
Z(jk± ~ 0.34. The bottom of the pocket is well fitted by a 
parabolic band ek± ~ (3ik\/2 + (32k2/2, where /3i and 02 
represent the inverse effective mass along and orthogonal 
to the nodal direction, respectively. At the highest dop- 
ing examined, x = 0.059, we obtain ^1 « 3.4, /32 ~ 1.8, 
which yields the anisotropy /3i//32 ~ 1-9. The value of 
the effective mass to* ~ l/v^AT^ converted to conven- 
tional units is TO* « 1.6TOe. This value is consistent with 
the mass measured in MQOJ^"— 



=0,;31=2.4,j62=2.2 



x=0.02,;ei=2.6,j52=2.2 




FIG. 6: (color online) Maps of the parity average disper- 
sion ~ (ek,+ + ek,-)/2 at different doping levels. Effec- 
tive mass along (/3i) and orthogonal (/32) to the nodal direc- 
tion is calculated by fitting the bottom of the pocket with 
Ek ~ fiik\/2 + /32k2/2. The ellipticity of hole pockets, repre- 
sented by y73i//32, increases with doping indicating a devia- 
tion from the rigid band approximation. 



Maps of the parity average dispersion, ek ~ (ck, 
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ek,-)/2, are shown in Fig. [6] for four values of doping. 
The total bandwidth « 2. 4 J « 0.3eV remains roughly 
the same independent of doping. This bandwidth is very 
close to the value indicated by ARPES in undoped sin- 
gle layer compound 4^— By comparing different doping 
levels in Fig. \6\ it is also clear that the ellipticity of the 
pocket increases with doping. More precisely, /3i is in- 
creasing while /32 is decreasing at higher doping. En- 
hancement of ellipticity with doping is consistent with a 
previous analytical calculation^^ The doping dependent 
ellipticity implies that the rigid band approximation is 
strictly speaking not valid. However, the change of el- 
lipticity, although being significant, is not dramatic, so 
the rigid band approximation is fairly reasonable. More- 
over, the average effective mass m* = 1/ •y//3i/32 is roughly 
doping independent. 

The present approach certainly supports the small 
Fermi surface Luttingcr's theorem. The doping calcu- 
lated according to Eq. ([TS)) coincides with the area of the 
small Fermi surface. It is worth noting that numerically 
this coincidence is quite nontrivial, the negative energy 
incoherent part of the Green's function is absolutely sig- 
nificant for this. 



IV. SPIN-CHARGE RECOMBINATION AND 
ARPES SPECTRAL FUNCTION 

In the SCBA approach spin and charge are separated, 
there are nonitinerant spins and there are itinerant spin- 
less holes. The hole has a pseudospin indicating a che- 
querboard magnetic sublattice, but the pseudospin is 
different from usual spin. In ARPES process spin and 
charge rccombine to physical electrons. To calculate 
the recombination amplitude we follow the approach of 
Ref. m modifying the approach to the bilayer case and 
to finite doping. ARPES measures electrons with true 
spin, regardless which sublattice the electron comes from. 
The annihilation operator of an electron from m-th plane 
(to = 1,2) is 



ie{t,]} 



— zk-r; 



(23) 



where a is the true spin. Notice that the definition ([23 
is properly normalized, because 



(0|4,i,^Ck,i,TlO) 



N 



2 

TV' 



(24) 



To establish the connection with hole operators in 
Eq. p^ . one needs to rotate electron operators into the 
fixed parity basis 



Ck,±,cr — (ck,l,(T ± Ck,2,( 



(25) 



The connection between electron and hole Green's func- 
tion is then associated with the following vertice o^^'^'^ 

flk = (0|d_k,+,tCk,+4|0) 



(0|/3q, + d-k-q, + ,4,Ck, + ,i|0) 



= (0|/3q. + d_k-q.-.iCk,-,i|0) 

&k,q,- = (0|/3q,_(i_k-q, + ,iCk,-,i|0) 

= (0|^q,_d_k-q,-,iCk, + ,i|0) 

Ck,q,-H = (OM_k+q, + ,iCk, + ,i"q,+ |0) 

= (0|d_k+q,-4Ck _4a]j_+|0) 

Ck,q- = (0|(i_k+q,+4Ck,-,ia]j _|0) 

= (0|(i_k+q_,iCk,-H,ia]j_|0) 



]^"q,-(26) 



Here |0) is the ground state of the doped system. The 
vertex Ok describes the process of instant creation of an 
electron with momentum k and a hole with momentum 
— k, |0) ^ cjj. ^ I'^^-k + fl*-*)' shown in Fig. [71(a). In 



(a) 



(b) 



(c) 



k -k-q k -k+q 



FIG. 7: Vertices involved in the photoemission process, as 
defined in Eq. ^ 



the figure electron is shown by the bold solid line mov- 
ing from the left, and the hole is shown by the double 
line. The direction of the electron line in Fig. [71 strictly 
speaking is not correct, as the electron has to appear in 
the final state only. Nevertheless, in figures we always 
show electron in the initial state just for a convenient 
graphical presentation. This way of presentation does 
not cause any problems since we always have only one 
electron. The vertex 6k, q, 7 describes the process of in- 
stant creation of an electron with momentum k, a hole 
with momentum — k — q and a magnon with momentum 
q. as shown in Fig. [7Kb). We define the parity index 7 
to be parity of the magnon. The vertex Ck,q.'y describes 
the process of annihilation of a magnon with momentum 
q and instant creation of an electron with momentum k 
and a hole with momentum —k+q, as shown in Fig.lTJc). 
Derivation of the vertices (j26p is presented in Appendix 
A. The analysis of ARPES in undoped parent compound 
at zero temperature requires only Ok and 5k, q vertices,— 
whereas the undoped compound at nonzero temperatures 
requires additional vertex Ck,qr^ because there are ther- 
mally excited magnons in the initial state. Presently we 



consider a doped compound at zero temperature, so there 
are no thermally excited magnetic fluctuations, but there 
are additional magnetic fluctuations due to doping. Thus 
all three vertices are involved in the present calculation. 

The electron Green's function is defined in the stan- 
dard way. 

/oo 
{Tcl±Jt)cK±AO)ydt (27) 
-oo 

Here a is the true spin index. Notice that the two spins 
are degenerate, so below we omit the spin index in Gc±- 
Dyson equations relating Gc and already calculated Gd 
are similar to that derived in Refs. [28II431 The equations 
are shown graphically in Fig. |S] and presented below in 
analytical form 



Ge+(e,k) 



aiGd+{e, -k) + k) + EL'i(e,k) 



.(1) 



2akGd+(e,-k) E^'+le^ k) + SL'i(e, k) 



(2) 



Gd+(e,-k) 



.(2) 



Gc-(e,k) 



(^,k) 

.(1) 



S]L'l(e,k) 

(1) 



aiGd- (e, -k) + S]Yi(e, k) + k) 
2akGd_ (e, -k) [s^'l (e, k) + (e, k) 



+ Gd_(e,-k) [l]V'l(e,k) + EL'|(e,k)] (28) 

We emphasize that the conventional self-energy T,^s ap- 
pears only in Dyson's equation (fTQ]) for Gd- Equations 
(PS)) contain different kinds of self-energy. The self-energy 

S^^y has two vertices ?'k,q.7 or Ck.q,7 denoted by the 

(2) 

crossed circle. The self-energy E^^' has only one vertex 
&k.q,7 or Ck,q.7 dcnotcd by the crossed circle and one ver- 
tex (7k, q (see Eq. (|15p ) shown by a simple attachment of 
the dashed line (magnon) to the double line (hole). The 
subscript {^6} in the self-energies shows parities of the 
intermediate magnon and hole, respectively, as displayed 
in Fig. ISl Expressions for Y]''^^ and T,''^g are presented in 



Appendix B. 

It is worth noting that meaning of Dyson's equations, 
Eq. (gSl) is different from that of Eq. ([15]). Eq. is a 
real dynamic equation for Gd±, which requires an itera- 
tive solution. On the other hand, Eq. (pS)) is just a rela- 
tion between Gc± and Gd± that describes the spin-charge 
recombination in the photoemission process. Once Gd± 
is known, we calculate the right hand side in Eq. pS]) and 
Gc± is determined. The Greens's function Gc± has the 
spectral representation similar to Eq. (jl7p and a spec- 
tral function similar to Eq. ([2^ . The spectral functions 
Sc± for X — 0.059 are displayed in Fig.[9l^a) and Fig.[9][b). 
The position of the quasiparticle peak is exactly the same 
as that in Gd±, as one expects from general considera- 
tions and it is also evident from Eqs. ([25)1 . However, the 
quasiparticle residue of Gc±, denoted by Zc±, is very dif- 
ferent from that of Gd±- Since Gc± represents electrons 
of either sublattice, there is no a Bloch theorem to require 
Zc± to be symmetric with respect to MBZ boundary. 



+ + + + + r ~\ / + 
= ^ — [®^^- — 

+ - 

+ - 

+ __- 




FIG. 8: (color online) Dyson's equations relating Gc± and 
Gd±- The thick line represents Gc±, and the double line rep- 
resents Gd±- 



In fact, our calculation shows a very asymmetric Z^^, 
as shown in Fig. [TO] along a certain line in BZ. Along 
the nodal direction, Zc± first increases from k = (0, 0) 
to the inner edge(red dotted vertical line) of the Fermi 
pocket. After that Zc± drops abruptly as k goes across 
the pocket inside the Fermi surface (although this part 
of the Green's function is not measurable by ARPES). 
A less steep drop of Zc± takes place from the outer edge 
(blue dotted vertical line) of the pocket to the point k = 
(7r,7r). The difference between Zc± on the inner side of 
the pocket and the outer side is very significant. 

The above analysis is performed within the double 
layer t - t' - t" - J model. The t - t' - t" - J model 
can originate from the single band Hubbard model or 
from a multi-band Hubbard model. The origin of the 
t — t' — t" — J model is not important for the dynamic 
equation ([TO]). However, the origin is important for the 
spin-charge recombination amplitude given by Eq. (|28p . 
Depending on the original model, there is an additional 
significant contribution to the asymmetry of the electron 
spectral function between inside and outside of MBZ. 
Analysis performed in Ref. [2^ shows that in the case of 
the single band Hubbard model the vertices in Eq. ([26]) 
should be modified by 



Ok 



Ok 



, J 
1 + -7. 
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FIG. 9: Spectral functions Sc+ (a) and Sc- (b) for different 
values of k along the nodal direction from (0,0) to (7r,7r). 
There is an offset 0.2 between subsequent values of k. The 
doping level is a; = 0.059. 



0.4 



0.3 



P: 0.2 



0.1 

















Gc-/ / 


i\ 












I I 
I 1 
1 1 


i\ 

i\ 
i\ 
i\ 






A\ 






1 1 

1 

1 1 


i\ 
iV 
1 






'/ '• 


i ' 
\ 
\ 






\ 


\ 

\\ 




n \ 
'/ • 
// ; 


■ 








\\ 
ft 











(0,0) 



(^,0) 



(;r/2,;r/2) {n,nl2) 



FIG. 10: (color online) The quasiparticle residue Sc along the 
same line as in Fig. O Doping level is a; = 0.059. Vertical 
dashed lines show Fermi points. 



J 



&k,q,± -5> ^k,q,± ( 1 + y7k 
Ck,q,± Ck,q,± ( 1 + Y7k 



(29) 



Following Ref. [2^ we call this correction the Hubbard 
model correction. Plots of the parity-averaged quasipar- 
ticle residue, [Z^^(k) + Z^(k)] /2, along the nodal di- 
rection are presented in Fig. [TT] for three different values 
of doping. In the same figure, we also present colour 
maps of the following function 



(52 



(31) 



(ek,+ - [if + (ek,- - [if 

This formula gives a way to image the quasiparticle 
residue at the Fermi surface. The broadening 8 can be 
taken result of disorder and/or experimental res- 
olution, or purely for the sake of imaging. We choose 
8 ~ 0.1, which corresponds to the halfwidth F = 2(5 = 
26mcV. 
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FIG. 11: (color online) Right column: The parity average 
residue [Z^(k) H- Z^(k)] /2 of the electron Green's func- 
tion along the nodal direction. Vertical dashed lines show 
Fermi points. The ARPES intensity is proportional to the 
residue outside the Fermi surface. Left column: colour maps 
of the ARPES intensity at the Fermi surface. The plots are 
presented for three values of doping, x = 0.02, 0.037, 0.059. 



which effectively modify quasiparticle residue Zc± by 



Ze± ^ Z^^ = Ze± ( 1 + -7k 



J 



(30) 



From the right column of Fig. [Tl] we observe that at 
the largest doping examined, x = 0.059, the ARPES in- 
tensity in the inner side (red dotted line) of the pocket 
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is about 4 times larger than that in the outer side (blue 
dotted line). Comparing plots in different doping levels, 
it is also clear that this asymmetry grows with doping. 
This asymmetry is also reflected in the left column of 
Fig. [TTJ where maps of Z ^ (k) at corresponding doping 
levels is presented. One clearly sees that the intensity of 
Z ^ (k) displays an arc shape along the Fermi surface. 

V. CONCLUSIONS 

At doping below 5.5% — 6% the bilayer cuprate 
YBa2Cu306+i/ is a collinear antiferromagnct. The 
doping-independent staggered magnetization at zero 
temperature is about 0.6/iB. This is the maximum value 
of magnetization allowed by quantum fluctuations of lo- 
calized spins. These experimental observations create a 
unique opportunity for theory to perform a controlled 
calculation of the electron spectral function at doping 
a; < 0.06. In the present work we perform such a calcu- 
lation within the framework of the extended t — J model 
with account of the Hubbard model corrections. The 
calculation employs the self-consistent Born approxima- 
tion (SCBA) that is parametrically justified because of 
the long range AF order with maximum possible stag- 
gered magnetization. To perform the work we have de- 
veloped/extended the SCBA to the finite doping case and 
to the bilayer system. The calculation clearly demon- 
strates that the Fermi surface consists of small hole pock- 
ets centered at (±7r/2, ±7r/2). This conclusion itself is 
a trivial one since we deal with the system with long 
range AF order with maximum possible staggered mag- 
netization. The small pocket is a direct consequence of 
Bloch theorem. What is nontrivial is that we quantify 
the asymmetry of the ARPES spectral function that is 
highly anisotropic at the Fermi surface. In particular at 
doping about 5% — 6% the ARPES intensity in the inner 
side of the pocket is about 4 times larger than that in 
the outer side. Overall the picture resembles Fermi arcs 
observed in ARPES. 

Our analysis shows that the hole band is not quite rigid 
under doping. In particular the ellipticity of the pocket 
increases with doping. However, the effect of changing 



ellipticity while being significant still is not dramatic, so 
the rigid band approximation is not that bad. The hole 
effective mass averaged over the Fermi surface is practi- 
cally doping independent, m* « l.Gnie. 

Our calculation demonstrates that due to a saddle 
point in the hole dispersion the electronic response of 
the system is peaked at q sa {n/2, 0), w « lOOmeV. These 
parameters are very close to those where the anomaly is 
observed in the breathing phonon mode. 

We also found that the antiferromagnetic correlations 
practically destroy the bilayer bonding-antibonding 
splitting. More precisely the correlations suppress the 
splitting by one order of magnitude. If without account 
of the correlations the splitting is about 200meV then at 
the doping 6% the value is reduced down to 20meV. 



VI. ACKNOWLEDGMENTS 

We thank P. Horsch, G. KhaliuUin, O. Jcpsen, A. I. 
Milstein, and A. Avella for stimulating discussions. A 
significant part of this work was done during stay of W.C. 
and O.P.S. at the Max Planck Institute for Solid State 
Research, Stuttgart, and stay of O.P.S. at the Yukawa 
Institute, Kyoto. W.C. and O.P.S. are very grateful to 
colleagues for hospitality and for stimulating atmosphere. 

The work was supported by ARC grant DP0881336. 
Computations for this project were performed using the 
National Computational Infrastructure, under project 
u66. 

This work was also supported by the Humboldt Foun- 
dation and by the Japan Society for Promotion of Sci- 
ence. 



Appendix A: calculation of hole-magnon vertices 

The hole-magnon vertex due to i?t,t^ in ^ contains 
two parts: the first part comes from in-plane nearest- 
neighbor hopping Ht- For example, the contribution of 
the 1st plane hopping in the case when all parities are 
positive is 



(0|aq,+<ik,+,t 
1 

iV(l/2 + m) 
t 

iV(l/2 + m) 



SlAUsl^ + (l-sl,]s-, 



|0) 



+ + 7k+ql'q.- 



(Al) 
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where we denote = + ^, and use the mean field decomposition cj ^cj ^Cj^i.fCij^^ — )• — cj j 4,Cj,i.t(cJ i ^Ci.i.f) = 
5^- ]^(l/2 + m). The second plane hopping iJj gives an equal contribution and altogether this results in the first term 

in gk,q,75 in Eq. (HH]). 

The contribution from the interlayer hopping is 

(0|aq,+dk,2,ti?t.4+c,i.i|0> 



-2ti 



iV(l/2 + m) 
iV(l/2 + m)' 



^(0|aq,+ 5^ e''^-. (Sr, + Sr^) |0) = t^y|(0|aq,+ {b^^^, + J |0) = 



|0) =*i\/]^(^iq,+ + ^'q,+ ) 



(A2) 



Rotating this to the parity basis, one recovers the second term in g^^q.-ys in Eq. (jlSp . 

The spin-charge recombination vertices (j26p are calculated in a similar way. The a-vertex is the following: 



Ok = (0|d_k,+,tCk,+4|0) 



(0| l^E ii.i^JM.i + E 424C».2,; j 



ivyT72 

The 6- vertex reads: 

''k,q,+ = (0|/3q, + (i_k-q. + ,iCk. + ,i|0) 

_ /o, /'q 



(A3) 



^yj72^ ^E4.r^'^''"'^^-^^^E^^M4e--'+E42,t^^^'^^''^-^^ E ^.24e-- I 10) 



/3q,H 



Similarly the c- vertex is: 



'q,+ 



(A4) 



Ck,q.+ = (0M-k+q, + 4Ck, + 4«l,+ |0) = y ^^^= « V ^^q^ + 



(A5) 



Appendix B: The hole self-energy 

Here we present explicit expressions for each self-energy in the Dyson's equations for Gd± and Gc± 



- y (^fek,q,,5k-q,q,.. - 0.^,, - X + ^0 ^ 7 ^'^''^ '^(^ + .q,, - - ,0 ' 

We remind that {7, 5} = {±, ±}. Notice that due to violation of cross-leg symmetry, the retarded and advanced part 
have different vertices, as shown in Fig.[51 see also Refs. [29ll30l It is worth mentioning that E-y^ is the true self-energy, 
while S^^y and T^i^j are dimensionless quantities describing the spin-charge recombination process. 
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Direct evaluation of (|Bip requires a three dimensional integration, two momenta and one frequency. This is too 
expensive computationally. Fortunately, using the Kramers-Kronig dispersion relation one can effectively remove one 
integration. Take S-i for example, its imaginary part is 

/mS;+_(e,k) 

^'^Jjo:^ -3k-q,q,+-y dxA^{x,k- q)S{e - ujq^+ - x) + gl_^^_ J dxB^{x,k - q}d{e + - x) 



(27r)2 



.9k-q.q. + -^- (e - '^q,+ . ^ - q)0{e - Wq 



(27r)2 L ^'--^■^ 
Imll^_ (e, k) + /mS^_ (e, k) 



5k,-q,+--B-(e + '^q,+ , k - q)6l(-e - Wq) 



(B2) 



Evaluation of ImT,^ and ImT,^ requires only a two-dimensional integration over momenta. The Kramers-Kronig 
relation is then applied to find the real part by the principal value integration 



i?eE^_(e,k) = - 



de'/mS]^_(e',k) r, r ^^ f de' Im^f_{e' ,-k) 



i?eSf_(e,k)=y ^ 



(B3) 
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